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Abstract.  An  electrophoretic  survey  of  42  populations  of  Anopheles  pseudopunctipennis  collected  throughout  its 
known  geographic  distribution  was  performed  to  clarify  the  taxonomic  status  of  this  important  malaria  vector  species. 
The  results  indicated  strong  differences  in  the  allele  frequencies  of  three  enzyme  loci  (glycerol  dehydrogenase,  6- 
phosphogluconate  dehydrogenase,  and  phosphoglucomutase)  of  the  33  loci  analyzed.  No  fixed  electromorphic  differ¬ 
ences  separate  the  populations  of  An.  pseudopunctipennis.  The  populations  of  An.  pseudopunctipennis  showed  little 
genetic  divergence,  with  Nei  distances  ranging  from  0  to  0.079.  A  comparison  of  An.  pseudopunctipennis  data  with 
either  one  of  three  other  Anopheles  species  showed  a  high  genetic  distance  of  0.335  with  a  closely  related  species. 
An.  franc iscanus’,  0.997  with  An.  crucians,  and  2.355  with  An.  (Nyssorhynchus)  albimanus.  Geographic  populations 
of  An.  pseudopunctipennis  were  classified  into  three  clusters;  one  cluster  included  populations  collected  in  North 
America  (United  States  and  Mexico)  and  Guatemala,  one  cluster  included  populations  from  Belize  and  South  America 
(Colombia,  Ecuador,  Peru,  Chile,  and  Argentina);  and  one  cluster  was  represented  by  populations  from  the  Island  of 
Grenada  (type-locality  of  An.  pseudopunctipennis).  Based  on  our  isozyme  analyses,  we  defined  these  clusters  as  three 
geographic  populations  of  An.  pseudopunctipennis.  Of  the  two  mainland  populations,  one  extends  from  the  southern 
United  States  south  through  Mexico  and  Guatemala,  and  the  other  extends  north  from  southern  South  America  through 
Central  America  to  Belize.  These  two  geographic  populations  converge  in  southern  Mexico  and  northern  Central 
America.  One  part  of  the  convergence  zone  was  identified  in  the  area  of  eastern  Guatemala  and  southern  Belize. 


Anopheles  pseudopunctipennis  Theobald  1901  is  a  major 
vector  of  human  malaria  in  the  foothills  and  mountainous 
areas  of  Mexico,  throughout  Central  America,  and  in  the 
Andean  countries  of  South  America.  It  is  often  the  only  vec¬ 
tor  present  in  areas  at  an  altitude  above  600  m.  It  is  found 
from  the  United  States  (south  of  40°N)  to  the  northern  part 
of  Argentina  (30°S)  along  the  Andes,  with  an  eastern  exten¬ 
sion  into  Venezuela  and  the  Lesser  Antilles.  In  the  most 
extensive  review  of  this  species,  Aitken'  asserted  that  “Be¬ 
cause  of  its  extensive  north  south  distribution,  its  inconsis¬ 
tencies  as  a  vector  of  malaria  in  the  Americas,  and  because 
of  the  conflicting  reports  concerning  its  morphology  and 
habits.  An.  pseudopunctipennis  has  come  to  be  an  extremely 
interesting  mosquito  to  study”.  Bruce-ChwatF  added  that  the 
extensive  geographic  distribution  and  the  high  level  of  vari¬ 
ability  of  this  species  have  led  to  speculation  that  An.  pseu¬ 
dopunctipennis  may  represent  a  complex  of  sibling  species. 
At  the  present  time,  the  taxonomic  status  of  An.  pseudo¬ 
punctipennis  is  extremely  important  because  of  its  wide  in¬ 
volvement  in  the  transmission  of  human  malarial  pathogens. 

Anopheles  pseudopunctipennis  was  first  described  by 
Theobald^  in  1901  from  Grenada  Island  (Lesser  Antilles). 
The  original  description  was  not  sufficient  for  accurate  iden¬ 
tification  of  the  species,  resulting  in  confusion  and  misiden- 
tifications  of  the  species.'  Between  1907  and  1912,  four  dif¬ 
ferent  names  were  applied  to  the  species  (three  are  currently 
in  synonymy).  From  1901  to  1950,  five  subspecies  and  one 


variant  of  An.  pseudopunctipennis  were  morphologically  de¬ 
scribed  from  different  areas  of  South  America.''  Anopheles 
franciscanus  was  a  synonym  of  An.  pseudopunctipennis  for 
28  years  and  then  considered  a  subspecies  of  An.  pseudo¬ 
punctipennis  for  another  40  years.  Only  in  1972  was  An. 
franciscanus  elevated  to  the  species  level. ^  In  1992  and 
1993,  Estrada-Franco  and  others^*  stated  that  An.  pseudo¬ 
punctipennis  constituted  a  complex  of  two  species.  An.  pseu¬ 
dopunctipennis  A,  a  species  from  central  Mexico,  and  An. 
pseudopunctipennis  B,  a  species  from  the  Andes  of  Peru  and 
Bolivia. 

Since  1991,  an  extensive  investigation  of  An.  pseudo¬ 
punctipennis  from  its  whole  geographic  range  has  been  un¬ 
dertaken  and  the  results  are  the  subject  of  this  report.  Our 
study  of  An.  pseudopunctipennis  was  a  genetic  analysis  by 
isozyme  electrophoresis  with  the  following  emphasis:  1)  an¬ 
alyzing  genetic  differentiation  and  variability  of  An.  pseu¬ 
dopunctipennis  within  its  known  geographic  range,  from  the 
type-locality  in  the  Caribbean  (Grenada  Island)  to  North, 
Central,  and  South  America;  2)  comparing  genetic  differ¬ 
entiation  and  variability  of  An.  pseudopunctipennis  popula¬ 
tions  that  correspond  to  the  five  subspecies  and  one  variant 
described  in  the  literature;  3)  comparing  genetic  differenti¬ 
ation  and  variability  of  An.  pseudopunctipennis  populations 
from  a  range  of  altitudes;  and  4)  comparing  the  genetic  pro¬ 
file  of  the  An.  pseudopunctipennis  populations  with  other 
associated  species  of  Anopheles. 
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Figure  1 .  Collecting  areas  of  Anopheles  pseudopunctipennis  in 
10  countries  throughout  the  species’  geographic  distribution.  Each 
dot  corresponds  to  a  collection  of  multiple  samples. 


MATERIALS  AND  METHODS 

Mosquito  populations.  Anopheles  pseudopunctipennis 
was  collected  as  wild  larvae  and  pupae  along  its  known  geo¬ 
graphic  range  (Figure  1)  and  at  altitudes  from  sea  level  up 
to  2,500  m.  The  countries  were  chosen  according  to  critical 
locations,  such  as  the  type-locality  (Grenada  Island),  and  lo¬ 
calities  where  the  five  different  subspecies  and  one  variant 
of  An.  pseudopunctipennis  were  described,"^  as  well  as  areas 
providing  a  spatial  representation  of  the  geographic  distri¬ 
bution  of  the  species.  In  South  America,  collections  were 
made  where  the  variant  bifoliata'^  occurs  in  Colombia,  where 
the  subspecies  levicastilloP°  and  rivadeneirai'^^'^  occur  in 
Ecuador,  where  the  subspecies  neghmei  and  noeP^  occur  in 
Chile,  and  where  the  subspecies  patersonp'^  occurs  in  Ar¬ 
gentina.  The  site-specific  collections  were  positive  for  all 
subspecies  and  variant  except  in  Chile  where  the  two  sub¬ 
species  of  An.  pseudopunctipennis  seem  to  have  been  erad¬ 
icated  from  the  two  locations  described  by  Mann.'^  However, 
we  found  populations  of  An.  pseudopunctipennis  near  Arica, 
in  the  northern  part  of  Chile,  approximately  200  km  from 
the  two  subspecies’  type-localities. 

All  An.  pseudopunctipennis  adults  tested  were  collected 
as  wild  larvae  or  pupae  and  reared  individually  to  the  adult 
stage.  Fourth  instar  larval  and/or  pupal  exuviae  were  pre¬ 
served  and  each  specimen  was  recorded.  An  adult  with  as¬ 
sociated  larval  and  pupal  skins  is  referred  to  as  a  taxonomic 
series.  Each  taxonomic  series  was  identified  to  species  and 
used  to  support  the  isozyme  study.  In  most  samples,  an  av¬ 
erage  of  40  adults  (50%  of  both  males  and  females)  were 
frozen  at  —  70°C  for  genetic  studies  and  some  adults  were 
pointed  on  pins  for  taxonomic  vouchers  and  deposited  at  the 
Smithsonian  Institution  (Washington,  DC). 

From  1991  to  1993,  we  collected  a  total  of  42  samples  in 
10  different  countries  represented  by  more  than  2,000  wild 
specimens  of  An.  pseudopunctipennis  (Table  1).  To  evaluate 
the  taxonomic  significance  of  observed  differences  among 
An.  pseudopunctipennis  populations,  three  other  Anopheles 


species  were  also  included  in  the  electrophoretic  data:  a 
closely  related  species.  An.  franciscanus  McCracken,  An. 
crucians  Wiedemann,  and  An.  (Nyssorhynchus)  albimanus 
Wiedemann.  Wild  specimens  of  An.  franciscanus  were  col¬ 
lected  near  Redding,  California,  and  both  An.  crucians  and 
An.  albimanus  were  collected  in  northern  Belize.  As  for  An. 
pseudopunctipennis,  all  three  Anopheles  species  were  col¬ 
lected  as  wild  larvae  and  pupae  and  reared  individually  to 
the  adult  stage.  An  average  of  40  adults  per  population  were 
stored  at  — 70°C  until  used  in  electrophoretic  studies. 

Electrophoresis.  Isozymes  were  separated  by  horizontal 
starch  gel  electrophoresis.  Staining  procedures  were  adapted 
from  the  reports  of  Harris  and  Hopkinson,'^  Selander  and 
others,'^  and  Shaw  and  Prasad.'"^  A  total  of  45  enzyme  sys¬ 
tems  were  screened  using  three  different  buffers:  the  lithium 
buffer  system'*  (LiOH,  pH  8.5),  the  morpholine  buffer  sys¬ 
tem'^  (morph,  pH  6.1),  and  the  Tris-citrate  buffer  system^® 
(TCss,  pH  6.7).  Of  the  45  enzyme  systems  tested,  25  showed 
good  allelic  resolution,  including  33  putative  loci  (Table  2). 
For  each  locus,  the  most  frequent  electromorph  was  desig¬ 
nated  the  100  allele  and  all  other  alleles  were  measured  rel¬ 
ative  to  it.  Additional  details  of  the  electrophoretic  procedure 
have  been  described  by  Manguin  and  others,^'  and  specific 
buffer  systems  and  staining  recipes  are  available  upon  re¬ 
quest  from  one  of  the  authors  (SM). 

Electromorph  genotype  frequencies  were  used  as  input  for 
the  computer  program  BIOSYS-l.^^  Analysis  of  each  pop¬ 
ulation  included  computation  of  allele  frequencies,  hetero¬ 
zygosity  per  locus,  additional  measures  of  genetic  variabil¬ 
ity,  and  a  test  for  conformance  with  the  Hardy- Weinberg 
equilibrium  at  single  loci  by  chi-square  analysis.  Differen¬ 
tiation  among  the  populations  was  measured  by  F-statistics. 
Nei’s^^  unbiased  and  Rogers’^''  genetic  distances  were  clus¬ 
tered  by  the  unweighted  pair  group  method  using  the  arith¬ 
metic  average  (UPGMA)  to  produce  the  phenogram. 

RESULTS 

A  total  of  42  samples,  representing  more  than  2,000  wild 
specimens  of  An.  pseudopunctipennis,  were  collected  in  10 
different  countries  of  the  Caribbean,  North,  Central,  and 
South  America  (Figure  1  and  Table  1).  Anopheles  pseudo¬ 
punctipennis  populations  within  each  country  showed  neg¬ 
ligible  genetic  differentiation  and  were  combined  to  form  a 
geographic  area.  A  total  of  12  geographic  areas  correspond¬ 
ing  to  each  of  the  10  countries,  including  northern  (Monter¬ 
rey)  and  southern  (Tapachula)  Mexico,  and  the  Pacific  to  the 
Atlantic  sides  of  Guatemala  were  compared  in  our  analyses. 

Heterozygosity.  The  isozyme  comparison  of  33  loci  (Ta¬ 
ble  2)  among  the  populations  of  An.  pseudopunctipennis 
from  North,  Central  and  South  America  indicated  that  the 
mean  heterozygosity  ranged  from  0.022  to  0.101,  with  an 
average  (±  SEM)  of  0.059  (±  0.020)  across  all  mainland 
populations  (Table  3).  Mean  heterozygosity  for  the  popula¬ 
tions  from  Grenada  Island  was  much  lower^^  with  a  value 
of  0.003,  which  is  in  accordance  with  the  usual  low  level  of 
genetic  variability  of  isolated  island  populations. 

Genetic  heterogeneity.  The  F-statistics  (Fgj),  a  measure 
of  the  amount  of  differentiation  among  subpopulations,^’ 
showed  an  average  value  of  0.375  and  a  mean  index  of  fix¬ 
ation  of  individuals  relative  to  the  total  of  subpopulations 
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Table  2 

Electrophoretically  detected  enzyme  systems  of  Anopheles  pseudopunctipennis 


Enzyme  system 

E.C.  number* 

Symbol 

No.  of  locit 

Buffert 

Aconitase 

4.2. 1.3 

AGON 

2 

TCss 

Adenylate  kinase 

2.7.4.3 

AK 

1 

TCss 

Aldehyde  oxidase 

L2.3.1 

AO 

1 

LiOH 

Arginine  kinase 

2.7.3.3 

ARGK 

1 

LiOH 

Esterase 

3. 1.1.1 

EST 

1 

Morph 

Fumarase 

4.2.L2 

FUM 

1 

TCss 

Glutamate  oxaloacetate  transaminase 

2.6. 1.1 

GOT 

2 

Morph 

Glutathione  reductase 

L6.4.2 

GR 

2 

TCss 

Glycerol  dehydrogenase 

1.1.1.72 

GCD 

1 

Morph 

a-glycerophosphate  dehydrogenase 

1. 1.1.8 

GPDH 

1 

Morph 

Glyceraldehyde-3-phosphate  dehydrogenase 

1.2.1.12 

G3PDH 

1 

Morph 

Hexokinase 

2.7. 1.1 

HK 

2 

Morph 

3-hydroxyacid  dehydrogenase 

1.1.1.30 

HAD 

1 

Morph 

Isocitrate  dehydrogenase 

1.1.1.42 

IDH 

2 

Morph 

Lactate  dehydrogenase 

1.1.1.27 

LDH 

1 

LiOH 

Leucine  amino  peptidase 

3.4.11.1 

LAP 

2 

LiOH 

Malate  dehydrogenase 

1.1.1.37 

MDH 

1 

Morph 

Malic  enzyme 

1.1.1.40 

ME 

1 

Morph 

Mannose-6-phosphate  isomerase 

5.3. 1.8 

MPI 

1 

Morph 

6-phosphogluconate  dehydrogenase 

1.1.1.44 

6PGD 

1 

Morph 

Phosphoglucomutase 

5.4.2.2 

PGM 

1 

Morph 

Phosphoglucose  isomerase 

5.3.L9 

PGI 

1 

LiOH 

Pyruvate  kinase 

2.7.1.40 

PK 

2 

TCss 

Sorbitol  dehydrogenase 

1.1.1.14 

SDH 

1 

LiOH 

Triose  phosphate  isomerase 

5.3.1. 1 

TPI 

2 

Morph 

*  Enzyme  commission  (E.C.)  number, 
t  Number  of  scorable  bands  per  phenotype. 

t  Refers  to  the  electrophoresis  buffer  (see  Materials  and  Methods). 


Table  3 

Measures  of  genetic  variation  of  Anopheles  pseudopunctipennis  (1-12)  and  An.  franciscanus  (13),  An.  crucians  (14),  and  An.  albimanus  (15)* 


Mean  heterozygosity 


Populations 

No.  of 
samples 

Mean  sample 
size/locus 

Mean  no.  of 
alleles/locus 

%  polymorphic 
locit 

Direct  count 

Hardy-Weinberg 

equilibrium 

(expected)+ 

1.  United  States 

1 

52.2 

1.6 

48.5 

0.101 

0.098 

(0.9) 

(0.1) 

(0.028) 

(0.027) 

2.  Monterrey,  Mexico 

2 

75.0 

2.2 

72.7 

0.066 

0.078 

(4.9) 

(0.2) 

(0.017) 

(0.020) 

3.  Tapachula,  Mexico 

4 

147.7 

2.5 

78.8 

0.061 

0.065 

(8.6) 

(0.2) 

(0.016) 

(0.017) 

4.  Pacific,  Guatemala 

4 

67.5 

2.0 

60.6 

0.044 

0.050 

(3.9) 

(0,2) 

(0.012) 

(0.014) 

5.  Atlantic,  Guatemala 

4 

191.2 

2.4 

72.7 

0.050 

0.052 

(5.4) 

(0.2) 

(0.013) 

(0.014) 

6.  Belize 

4 

105.7 

1.8 

66.7 

0.074 

0.069 

(6.8) 

(0.1) 

(0.032) 

(0.025) 

7.  Grenada 

2 

69.5 

1.1 

12.1 

0.003 

0.003 

(3.6) 

(0.1) 

(0.002) 

(0.002) 

8.  Colombia 

1 

20.3 

1.5 

45.5 

0.059 

0.067 

(0.6) 

(0.1) 

(0.017) 

(0.021) 

9.  Ecuador 

4 

87.7 

1.7 

51.5 

0.069 

0.072 

(3.6) 

(0.1) 

(0.023) 

(0.024) 

10.  Peru 

5 

125.8 

1.9 

57.6 

0.039 

0.043 

(7.3) 

(0.2) 

(0.011) 

(0.013) 

11.  Chile 

6 

96.6 

1.5 

39.4 

0.022 

0.021 

(4.6) 

(0.1) 

(0.008) 

(0.008) 

12.  Argentina 

5 

80.1 

1.8 

51.5 

0.061 

0.065 

(4.9) 

(0.2) 

(0.021) 

(0.023) 

13.  An.  franciscanus 

1 

9.5 

1.4 

33.3 

0.084 

0.088 

(0.3) 

(0.1) 

(0.027) 

(0.028) 

14.  An.  crucians 

1 

32.2 

1.5 

39.4 

0.078 

0.071 

(2.2) 

(0.1) 

(0.029) 

(0.025) 

15.  An.  albimanus 

1 

35.8 

1.6 

30.3 

0.051 

0.057 

(2.5) 

(0.2) 

(0.023) 

(0.026) 

*  Values  in  parentheses  are  standard  errors. 

t  A  locus  is  considered  polymorphic  if  the  frequency  of  the  most  common  allele  does  not  exceed  0.99. 
t  Unbiased  estimate. 
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Table  4 

F-statistics  analysis  of  polymorphic  loci  in  42  populations  of  Anoph¬ 
eles  pseudopunctipennis* 


Locus 

F,s 

Fit 

Fst 

Acon-l 

0.020 

0.047 

0.027 

Acon-2 

-0.006 

-0.002 

0.004 

Ak-3 

0.016 

0.056 

0.041 

Ao 

-0.019 

-0.002 

0.017 

Argk 

0.330 

0.348 

0.027 

Est 

0.140 

0.363 

0.259 

Fum 

0.075 

0.089 

0.015 

Gcd 

0.084 

0.776 

0.755 

Got-] 

-0.031 

-0.007 

0.023 

Got-2 

0.017 

0.152 

0.138 

Gpdh 

-0.071 

-0.010 

0.057 

G3pdh 

-0.000 

0.019 

0.019 

Gr 

0.051 

0.231 

0.189 

Had 

0.008 

0.053 

0.045 

Hk-1 

0.212 

0.390 

0.226 

Hk-2 

0.015 

0.032 

0.017 

Idh-1 

-0.383 

-0.057 

0.235 

ldh-2 

0.072 

0.146 

0.079 

Lap-1 

0.082 

0.133 

0.055 

Lap-2 

0.037 

0.068 

0.032 

Ldh 

-0.151 

-0.019 

0.115 

Mdh 

-0.009 

-0.003 

0.006 

Me 

0.169 

0.176 

0.009 

Mpi 

0.205 

0.255 

0.063 

6Pgd 

0.100 

0.618 

0.576 

Pgi 

-0.004 

-0.001 

0.003 

Pgm 

0.069 

0.635 

0.608 

Pk-1 

-0.016 

-0.006 

0.011 

Pk-2 

-0.050 

0.027 

0.073 

Sdh 

-0.053 

0.111 

0.156 

Tpi-1 

-0.075 

-0.007 

0.063 

Tpi-2 

0.178 

0.204 

0.031 

Mean 

0.036 

0.397 

0.375 

*F|s  =  fixation  indices  of  individuals  relative  to  the  total  subpopulations;  Fn-  =  fixation 
indices  of  individuals  relative  to  the  total  populations;  F^t  =  F  statistics.  For  definitions  of 
loci,  see  Table  2. 


(F,s)  value  of  0.036  when  all  An.  pseudopunctipennis  pop¬ 
ulations  were  analyzed  (Table  4).  Three  loci,  glycerol  de¬ 
hydrogenase  {Gcd),  6-phosphogluconate  dehydrogenase 
{6Pgd),  and  phosphoglucomutase  {Pgm),  among  a  total  of 
33  showed  great  differentiation,  with  values  of  0.755,  0.576, 
and  0.608,  respectively.  In  the  case  of  Gcd  (Figure  2A),  the 
populations  of  An.  pseudopunctipennis  from  North  America 
(United  States  and  Mexico)  and  Guatemala  have  a  high  fre¬ 
quency  (92%)  for  allele  Gcd^22  (Table  5),  whereas  all  the 
other  An.  pseudopunctipennis  populations  showed  a  high  fre¬ 
quency  (90-100%)  for  allele  Gc<i,oo-  For  6Pgd  (Figure  2B), 
all  An.  pseudopunctipennis  populations  have  a  high  frequen¬ 
cy  (84-98%)  for  allele  dPgd^^Q,  except  the  ones  from  Gre¬ 
nada,  which  have  a  very  high  frequency  (100%)  for  allele 
6Pgd]2%-  With  Pgm  (Figure  2C),  North  America,  Guatemala, 
and  Grenada  have  a  high  frequency  (90-100%)  for  allele 
Pgm^QQ,  whereas  all  the  An.  pseudopunctipennis  populations 
from  South  America  and  Belize  have  a  high  frequency 
(87%)  for  allele  Pgm, 23.  No  fixed  differences  for  the  three 
loci  {Gcd,  6Pgd,  and  Pgm)  or  any  other  loci  have  been  found 
for  all  An.  pseudopunctipennis  populations  (Table  5). 

Electromorph  frequency  data  for  the  33  enzyme  loci  stud¬ 
ied  with  all  An.  pseudopunctipennis  populations  are  shown 
in  Appendix  A.  Only  one  locus,  glutathione  reductase-2,  was 
monomorphic  in  all  different  populations  and  species.  No 


electrophoretic  activity  was  shown  by  An.  crucians  for  iso¬ 
citrate  dehydrogenase-2  and  by  An.  albimanus  for  leucine 
amino  peptidase-2,  mannose-6-phosphate  isomerase  {Mpi), 
and  sorbitol  dehydrogenase.  Significant  departures  from  Har¬ 
dy- Weinberg  expectations  were  observed  in  only  10  cases 
(loci  indicated  by  $  in  Appendix  A)  among  the  396  com¬ 
parisons  (2.5%).  The  resolution  of  esterase  and  Mpi  was 
sometimes  poor,  which  might  have  introduced  some  scoring 
errors.  Five  populations  demonstrated  a  deficiency  of  het¬ 
erozygotes  from  expected  proportions  in  either  arginine  ki¬ 
nase,  fumarase,  Gcd,  glyceraldehyde-3-phosphate  dehydro¬ 
genase,  malic  enzyme,  Mpi,  or  Pgm.  Since  rare  alleles  were 
involved,  little  meaning  can  be  attached  to  these  significant 
(probability  <  1%)  deviations. 

Genetic  structure  of  subpopulations.  From  the  frequen¬ 
cy  data,  both  Nei’s'^  unbiased  and  Rogers’'**  distance  matri¬ 
ces  were  calculated  for  the  different  populations  and  species 
(Table  6).  Nei’s'^  unbiased  distance  was  chosen  for  ease  in 
comparing  these  results  with  those  of  Estrada-Franco  and 
others’  and  were  clustered  by  UPGMA  to  produce  the  phen- 
ograms  shown  in  Figures  3  and  4.  Phenograms  produced 
using  other  distance  measures,  such  as  Nei’s,’*  Rogers,’** 
modified  Rogers,”  Cavalli-Sforza  and  Edwards  chord  and 
arc,’^  produced  nearly  identical  branching  patterns.  The  three 
major  groupings  shown  on  Eigure  3  were  produced  by  each 
of  the  methods  listed  above. 

Populations  of  three  An.  pseudopunctipennis  subspecies, 
levicastilloi  and  rivadeneirai  from  Ecuador,  and  patersoni 
from  Argentina,  and  the  variant  bifoliata  from  Colombia, 
were  electrophoretically  compared  with  the  other  popula¬ 
tions  from  South  America.  The  results  showed  no  significant 
differences  among  the  populations. 

The  comparison  of  all  populations  of  An.  pseudopuncti¬ 
pennis  showed  some  differences  in  the  allele  frequencies, 
but  the  Nei’s  index  of  genetic  distance  indicated  a  high  de¬ 
gree  of  similarity,  with  values  ranging  from  0  to  0.079  (Table 
6).  The  phenogram  showed  three  clusters  of  An.  pseudo¬ 
punctipennis  populations  (Figure  3):  one  cluster  from  North 
America  (United  States  and  Mexico)  and  Guatemala,  a  sec¬ 
ond  cluster  from  Belize  and  South  America  (Colombia,  Ec¬ 
uador,  Chile,  Peru,  and  Argentina),  and  a  third  cluster  rep¬ 
resented  by  populations  from  Grenada  only.  The  phenogram 
indicated  that  An.  pseudopunctipennis  populations  from 
North  America  and  Guatemala  formed  a  group  with  very 
low  genetic  distance  (less  than  0.010).  The  F^t  analysis  using 
Wright’s  categories”  of  the  subpopulations  represented  by 
the  North  America-Guatemala  group  indicated  a  negligible 
differentiation,  with  a  mean  value  of  0.049.  The  second 
group  included  populations  from  Belize  and  South  America 
with  genetic  distances  less  than  0.024.  The  F^j  analysis  for 
populations  from  South  America  and  from  Belize-South 
America  showed  moderate  differentiation,  with  mean  values 
of  0.186  and  0.192,  respectively.  The  third  group,  repre¬ 
sented  by  populations  from  Grenada,  had  a  genetic  distance 
at  a  level  of  0.060  if  compared  with  populations  from  North, 
Central,  and  South  America.  The  two  major  differences  be¬ 
tween  the  populations  of  Grenada  and  all  the  others  were  1) 
the  different  allele  frequencies  for  6Pgd  and  2)  the  lack  of 
heterozygotes  due  to  isolation  of  the  Grenada  populations. 

Genetic  structure  compared  with  other  Anopheles  spe¬ 
cies.  The  comparison  of  An.  pseudopunctipennis  populations 
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Figure  2.  Electrophoretic  pattern  of  three  enzyme  loci  in  Anopheles  pseudopunctipennis  populations  from  Mexico,  Peru,  Belize,  and 
Grenada.  A,  Gcd  =  glycerol  dehydrogenase.  B,  6Pgd  =  6-phosphogluconate  dehydrogenase.  C,  Pgm  =  phosphoglucomutase.  Values  on  the 
right  are  the  relative  mobilities  of  alleles. 


with  either  of  the  three  Anopheles  species.  An.  franciscanus. 
An.  crucians,  and  An.  albimanus,  showed  increasingly  higher 
genetic  distances  (Figure  4).  When  the  comparison  involves 
two  species  of  the  subgenus  Anopheles,  the  genetic  distance 
varies  from  0.335  with  An.  franciscanus,  a  closely  related 
species  of  An.  pseudopunctipennis,  to  0.997  with  An.  cru¬ 
cians.  In  the  case  of  the  comparison  with  An.  albimanus, 
which  belongs  to  the  Nyssorhynchus  subgenus,  the  genetic 
distance  is  much  higher,  with  a  value  of  2.355. 

These  comparisons  emphasize  the  great  similarity  existing 
among  the  populations  of  An.  pseudopunctipennis  from  all 
geographic  areas. 

DISCUSSION 

Results  of  this  study  on  An.  pseudopunctipennis  indicate 
that  strong  allele  frequency  differences  occurred  for  three  of 


the  33  loci  {Gcd,  6Pgd,  and  Pgm),  but  no  fixed  differences 
were  found  for  any  of  these  loci.  The  mean  heterozygosity, 
which  reflects  the  genetic  variability,  was  very  low  for  An. 
pseudopunctipennis  populations  on  Grenada  Island  with  a 
value  of  0.003  (Table  3).  In  contrast,  the  mean  heterozygos¬ 
ity  of  mainland  An.  pseudopunctipennis  populations  varied 
from  0.022  to  0.101,  with  an  average  (±  SEM)  of  0.059  (± 
0.020).  The  mean  ±  SEM  heterozygosity  for  An.  pseudo¬ 
punctipennis  seemed  to  be  lower  than  that  expected  (0.115 
±  0.009)  for  other  Diptera  groups. 

In  the  case  of  An.  pseudopunctipennis  populations  from 
Grenada  and  Chile,  the  lack  of  alleles  for  6Pgd  is  part  of 
the  general  phenomenon  of  low  allelic  polymorphisms  found 
in  these  two  populations  (Table  3).  These  findings  are  com¬ 
patible  with  the  general  observation  that  species  or  popula¬ 
tions  that  are  distributed  over  a  variety  of  environmental 
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Table  5 

Relative  allele  frequencies  for  three  loci  of  Anopheles  pseudopunctipennis  clustered  in  three  geographic  groups  and  An.  franciscanus,  An. 
crucians,  and  An.  albimanus* 


Locust 

Gcd 


6Pgd 


Pgm 


An.  pseudopunctipennis 


Allelet 

United  States, 
Mexico, 
Guatemala 

Belize, 

South  America 

Grenada 

An.  franciscanus 

An.  crucians 

An.  albimai 

United  States 

Belize 

Belize 

n 

271 

299 

32 

10 

4 

3 

250 

0 

0 

0 

0 

0 

1.000 

173 

0.004 

0 

0 

0 

1.000 

0 

142 

0.020 

0 

0 

0 

0 

0 

122 

0.917 

0.003 

0 

0 

0 

0 

100 

0.057 

0.901 

1.000 

0.950 

0 

0 

77 

0.002 

0.095 

0 

0 

0 

0 

53 

0 

0 

0 

0.050 

0 

0 

(H) 

0.148 

0.144 

0 

0.100 

0 

0 

n 

617 

618 

95 

10 

39 

43 

225 

0 

0 

0 

0 

0 

0.012 

191 

0 

0 

0 

0 

0 

0.965 

147 

0.008 

0.014 

0 

0 

0 

0.023 

128 

0.002 

0.014 

1.000 

0 

0 

0 

100 

0.981 

0.839 

0 

1.000 

0.974 

0 

54 

0.006 

0.121 

0 

0 

0.026 

0 

23 

0.004 

0.013 

0 

0 

0 

0 

(H) 

0.039 

0.194 

0 

0 

0.051 

0.070 

n 

616 

609 

95 

10 

40 

42 

151 

0.006 

0.002 

0 

0 

0 

0.976 

123 

0.075 

0.865 

0 

0 

0 

0.012 

100 

0.899 

0.126 

1.000 

0.900 

0.975 

0.012 

75 

0.020 

0.006 

0 

0.100 

0.025 

0 

(H) 

0.172 

0.163 

0 

0.200 

0.050 

0.048 

*Bold  numbers  indicate  the  highest  frequency  of  each  locus  and  each  population, 

+Gcii  =  glycerol  dehydrogenase;  6Pgd  =  6-phosphogluconate  dehydrogenase;  Pgm  =  phosphoglucomutase, 
tn  =  no.  of  specimens;  H  =  heterozygosity  (direct  count)  per  locus. 


conditions  are  most  likely  to  be  genetically  heterozygous  or 
polymorphic.  Conversely,  species  or  populations  limited  in 
distribution  or  restricted  to  special  habitats  are  less  poly¬ 
morphic.^®  The  Grenada  Island  populations  differed  from 
mainland  populations  by  having  allele  frequency  differences 
for  6Pgd  and  by  lacking  heterozygotes  due  to  geographic 
isolation.  The  Grenada  Island  populations  showed  a  low  al¬ 
lelic  polymorphism  of  12.1%  compared  with  the  mainland 
populations,  with  values  ranging  from  39.4%  to  78.8%.  Av¬ 
erage  heterozygosity  over  the  33  loci  examined  was  also  low 
(3%),  reflecting  founder  effects.^'  The  relatively  higher  Nei 
genetic  distance  obtained  for  the  Grenada  populations 
(0.079)  compared  with  mainland  populations  (Table  6)  is 


partly  due  to  differences  in  allelic  polymorphism.  In  partic¬ 
ular,  the  absence  of  rare  alleles  is  an  indication  that  a  major 
reduction  in  the  gene  pool  occurred  in  the  recent  evolution¬ 
ary  history  of  this  island  population. This  phenomenon  can 
result  from  a  drastic  reduction  in  population  size  followed 
by  population  expansion  from  a  small  number  of  individuals. 
Such  fluctuations  in  population  size  are  common  occurrences 
in  insect  colonies. The  significant  loss  of  alleles  and  geo¬ 
graphic  isolation  of  An.  pseudopunctipennis  populations  of 
Grenada  are  limiting  factors  for  gene  flow  between  popula¬ 
tions  of  the  island  and  the  mainland.  As  a  result.  An.  pseu¬ 
dopunctipennis  populations  of  Grenada  may  eventually  form 
a  distinct  species  through  allopatric  speciation.^^  However, 


Table  6 

Matrix  of  genetic  distance  of  Anopheles  pseudopunctipennis* 


Population 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

1. 

United  States 

_ 

0.051 

0.047 

0.051 

0.050 

0.127 

0.121 

0.109 

0.102 

0.113 

0.112 

0.124 

2. 

Monterrey,  Mexico 

0.008 

- 

0.028 

0.027 

0.032 

0.098 

0.090 

0.089 

0.077 

0.083 

0.082 

0.102 

3. 

Tapachula,  Mexico 

0.008 

0.001 

- 

0.017 

0.017 

0.094 

0.089 

0.097 

0.078 

0.087 

0.081 

0.105 

4. 

Pacific,  Guatemala 

0.009 

0.001 

0.000 

- 

0.014 

0.091 

0.082 

0.098 

0.077 

0.085 

0.077 

0.103 

5. 

Atlantic,  Guatemala 

0.009 

0.002 

0.000 

0.000 

- 

0.092 

0.082 

0.095 

0.074 

0.084 

0.081 

0.102 

6. 

Belize 

0.067 

0.049 

0.053 

0.054 

0.056 

- 

0.093 

0.075 

0.067 

0.051 

0.041 

0.069 

7. 

Grenada 

0.079 

0.057 

0.062 

0.061 

0.063 

0.064 

- 

0.087 

0.078 

0.073 

0.070 

0.076 

8. 

Colombia 

0.052 

0.040 

0.049 

0.055 

0.054 

0.027 

0.059 

- 

0.062 

0.055 

0.058 

0.066 

9. 

Ecuador 

0.045 

0.029 

0.034 

0.035 

0.036 

0.017 

0.046 

0.022 

- 

0.054 

0.053 

0.068 

10. 

Peru 

0.067 

0.049 

0.055 

0.058 

0.060 

0.010 

0.057 

0.020 

0.014 

- 

0.025 

0.030 

11. 

Chile 

0.067 

0.048 

0.054 

0.056 

0.059 

0.009 

0.062 

0.021 

0.014 

0.002 

- 

0.042 

12. 

Argentina 

0.076 

0.059 

0.067 

0.069 

0.072 

0.018 

0.052 

0.026 

0.021 

0.004 

0.011 

- 

*Values  above  diagonal  are  Rogers^'*  genetic  distances;  values  below  diagonal  are  Nei^’  unbiased  (1978)  genetic  distances. 
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Figure  3.  Unweighted  pair  group  method  using  the  arithmetic 
average  phenogram  from  Nei’s'^  unbiased  genetic  distance  matrix 
for  all  Anopheles  pseudopunctipennis  populations  (cophenetic  cor¬ 
relation  =  0.936). 
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Figure  4.  Unweighted  pair  group  method  using  the  arithmetic 
average  phenogram  from  Nei’s'^  unbiased  genetic  distance  matrix 
for  all  Anopheles  pseudopunctipennis  populations  and  An.  francis¬ 
canus,  An.  crucians,  and  An.  albimanus  (cophenetic  correlation  = 
0.990). 


the  decreasing  presence  of  An.  pseudopunctipennis  on  the 
island  of  Grenada,  possibly  due  to  reduced  genetic  variabil¬ 
ity  accompanied  by  a  diminished  capacity  of  the  species  to 
adapt  to  environmental  changes, might  result  in  the  species 
disappearing  from  Grenada. 

The  collection  site  in  northern  Chile  also  represented  a 
restricted  ecosystem.  The  Tarapaca  region  where  An.  pseu¬ 
dopunctipennis  was  collected  in  Chile  is  bordered  by  the 
Pacific  Ocean  to  the  west,  the  Andes  to  the  east,  and  deserts 
to  the  north  and  south.  Anopheles  pseudopunctipennis  was 
collected  along  two  rivers  flowing  from  the  Andes,  the  Rio 
Azapa  and  the  Rio  Lluta,  and  was  much  more  abundant 
along  the  latter  river.  Because  of  geographic  barriers.  An. 
pseudopunctipennis  populations  in  Chile  had  more  restricted 
distributions  and  possibly  more  reduced  gene  flow  than  any 
other  populations  on  the  continent.  In  such  a  restricted  en¬ 
vironment,  the  bottleneck  effect  is  expected  to  produce  pop¬ 
ulations  less  polymorphic. 

Insect  colonies  are  characterized  by  rapid  genetic  changes, 
such  as  loss  of  heterozygosity. Loss  of  population  hetero¬ 
zygosity  in  colonies  probably  reflects  what  occurs  in  restrict¬ 
ed  and  isolated  environments  with  the  gradual  loss  of  poly¬ 
morphic  alleles.  Although  rare  alleles  contribute  little  to  a 
natural  population’s  level  of  heterozygosity,  these  are  the 
very  alleles  that  are  likely  to  be  missing  from  a  colony.^^ 
The  decrease  or  absence  of  rare  alleles  is  a  small  but  sen¬ 
sitive  indicator  of  the  more  general  phenomenon  of  loss  in 
genetic  heterogeneity,  particularly  in  species  with  high  ge¬ 
netic  variability.^^ 

Overall  genetic  distance  within  An.  pseudopunctipennis 
populations  was  low,  with  values  ranging  from  0  to  0.079. 
These  genetic  distances  are  markedly  lower  than  the  value 
of  0.16  suggested  by  Avise^^  as  the  lower  limit  for  conspe- 
cific  populations.  Values  of  Fgj  (Table  4)  show  great  differ¬ 
entiation  among  An.  pseudopunctipennis  populations.  This 
level  of  differentiation  is  due  to  low  heterozygosities  of  pop¬ 
ulations  from  two  specific  countries,  Grenada  and  Chile;  dif¬ 
ferences  in  the  allele  frequency  of  three  loci  among  33;  and 
the  broad  geographic  spread  of  sample  sites.  Regardless,  the 
moderate  mean  of  Fjs  with  a  value  of  0.036  suggests  that 


random  mating  among  the  populations  of  An.  pseudopunc¬ 
tipennis  is  occurring^’  (Table  4), 

As  clearly  indicated  in  the  phenogram  (Figure  3),  An. 
pseudopunctipennis  populations  are  clustered  into  three  ge¬ 
netically  associated  groups:  1)  North  America  and  Guate¬ 
mala,  2)  Belize  and  South  America,  and  3)  Grenada  Island. 
A  portion  of  the  convergence  zone  or  suture  zone^®  where 
geographic  populations  1  and  2  meet  is  located  in  the  300- 
km  area  around  the  southern  border  between  Guatemala  and 
Belize  (Figure  5).  Where  the  convergence  zone  is  located 
outside  of  this  area  is  not  known.  The  interface  zone  be¬ 
tween  the  two  populations  may  be  static  and  define  where 
the  North  and  South  American  populations  came  together 
after  a  long  period  of  isolation.  Alternatively,  the  conver¬ 
gence  zone  may  represent  only  the  most  recent  location  of 
a  mobile  border  between  two  merging  geographic  popula¬ 
tions. 

Our  findings  have  some  similarities  with  the  results  of 
Estrada-Franco  and  others,^  who  found  “clear  distinctions 
between  populations  from  South  America  and  Mexico  at  two 
loci,  Gcd  and  Pgm"\  However,  unlike  their  statement  that 
“the  major  contribution  to  genetic  divergence  is  the  result 
of  fixed  differences  in  the  Gcd  and  Pgm  loci  between  Mex¬ 
ico  and  South  America,”  we  found  no  fixed  differences  ei¬ 
ther  at  these  or  other  loci.  Comparisons  of  Nei’s  genetic 
distances  show  that  between  populations  of  An.  pseudo¬ 
punctipennis  and  An.  franciscanus,  values  are  at  the  level  of 
0.335,  and  among  An.  pseudopunctipennis  populations,  val¬ 
ues  are  less  than  0.08.  The  differences  in  the  genetic  dis¬ 
tances  between  An.  pseudopunctipennis  and  An.  franciscan¬ 
us  populations  confirms  the  separation  of  the  two  closely 
related  species  and  reinforces  our  finding  of  genetic  homo¬ 
geneity  among  An.  pseudopunctipennis  populations. 

Our  study  represents  an  advance  over  earlier  work  on  pop¬ 
ulation  genetics  of  An.  pseudopunctipennis.^^  We  used  25 
enzyme  systems  for  which  were  scored  33  putative  loci  that 
were  used  to  compare  An.  pseudopunctipennis  populations. 
A  large  number  of  sites  were  sampled,  providing  42  samples 
spread  from  the  northern  to  southern  limits  of  the  species’ 
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Figure  5.  Convergence  zone  between  Anopheles  pseudopunctipennis  populations  from  Guatemala  and  Belize. 
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distribution,  with  an  eastern  extension  in  the  Caribbean. 
Samples  were  obtained  from  one  Caribbean,  two  North 
American,  two  Central  American,  and  five  South  American 
countries.  Populations  from  Grenada,  the  type-locality  of  An. 
pseudopunctipennis  species,  were  compared  with  mainland 
populations.  This  was  an  important  part  of  the  study  since 
specimens  from  Grenada  represent  the  true  An.  pseudopunc¬ 
tipennis.  Additionally,  heterozygosity  levels  were  measured 
for  An.  pseudopunctipennis  populations  as  well  as  for  three 
other  Anopheles  species. 

In  conclusion,  our  study  showed  distinct  differences  in  the 
allele  frequencies  of  Gcd  and  Pgm  between  North  and  South 
America,  although  we  distinguished  overlapping  frequen¬ 
cies.  A  comparison  of  An.  pseudopunctipennis  populations 
with  three  other  Anopheles  species.  An.  franciscanus.  An. 
crucians  (both  of  the  subgenus  Anopheles),  and  An.  albi- 
manus  (subgenus  Nyssorhynchus),  provided  a  perspective  for 
interpretation  of  Nei’s  genetic  distances. 

Based  on  the  evidence  of  our  isozyme  analyses,  the  42 
samples  of  An.  pseudopunctipennis  were  clustered  into  three 
geographic  populations  represented  by  1)  North  America 
(United  States  and  Mexico)  and  Guatemala,  2)  Belize  and 
South  America  (Colombia,  Ecuador,  Peru,  Chile,  and  Argen¬ 
tina),  and  3)  Grenada  Island.  Of  the  two  mainland  popula¬ 
tions  (1  and  2),  one  extends  from  the  southern  United  States 
through  Mexico  and  Guatemala,  and  the  other  extends  north 
from  southern  South  America  through  Central  America  to 
Belize.  These  two  geographic  populations  converge  in  south¬ 
ern  Mexico  and  northern  Central  America.  Part  of  a  con¬ 
vergence  zone,  situated  at  the  vicinity  of  the  southern  border 
between  Belize  and  Guatemala,  was  defined  for  the  two 
mainland  geographic  populations  of  An.  pseudopunctipennis. 
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t  n  =  no.  of  specimens;  H  =  heterozygosity  (direct  count)  per  locus;  Negative  values  indicate  cathodally  migrating  alleles, 
t  Locus  deviating  from  Hardy-Weinberg  equilibrium. 


